Absence of singular superconducting fluctuation corrections to thermal conductivity 



(N 
O 

o 



O 



-a 



Douglas R Niven and Robert A Smith 
School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, ENGLAND 



We evaluate the superconducting fluctuation corrections to thermal conductivity in the nor- 
mal state which diverge as T approaches T c . We find zero total contribution for one, two and 
three-dimensional superconductors for arbitrary impurity concentration. The method used is di- 
agrammatic many-body theory, and all contributions - Aslamazov-Larkin (AL), Maki- Thompson 
(MT), and density-of-states (DOS) - are considered. The AL contribution is convergent, whilst the 
divergences of the DOS and MT diagrams exactly cancel. 
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^vq ' The discovery of the high-T c superconductors has led to a renewed interestEJ in superconducting fluctuation 

corrections to normal statjj-transport propertiesB. Whilst much of the work has focused on the electrical resistivity, 
~ 1, p, several experimentsEfuu'El hamjpnaijtfid fluctuation corrections to the thermal conductivity, k. Since there is 
some dispute between thcoristsQ&BliljOO as to the predicted magnitude of the effect, we have performed a detailed 



some dispute between theoristsDH'oEirOiij as to the predicted magnitude of the effect, we have performed a detailed 
£j . microscopic calculation valid for all impurity concentrations. We find no divergent fluctuation contribution, and 
^ 1 conclude that the experimental features seen near T c must have some other physical origin. 

Let us try to understand the reason for the lack of singular fluctuation contributions to thermal conductivity. 
There are several processes involved, and we will try to develop a physical pictureEj for each. The Aslamazov-Larkin 
| (AL) process involves the transfer of heat by fluctuation Cooper pairs. The corresponding term for the electrical 
+^ ■ conductivity has the strong divergence 

oal ~ (T - T c ) d ' 2 - 2 . (1) 



The size of the contribution to thermal conductivity can be estimated from Eq. ([!]) using the Wiedemann-Franz law, 
q . which has the general form 

o : 2 

KT-fMM a, (2) 
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00 • where UbTq is the amount of heat, and Qo the electric charge, carried by the excitations in a given system. For 
£f~) ' fluctuation Cooper pairs, T ~ T — T c , and Q = 2e so that 

o 

© [ ~ (^r^) — ~ (t - T c yi\ ( 3) 

which is clearly non-singular as T — > T c , The density-of-states (DOS) correction arises from the fact that when 
electrons form fluctuation Cooper pairs, they cannot simultaneously act as normal electrons; there is a corresponding 
ei ' decrease in the normal state density-of-states and hence normal state thermal conductivity 
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O ■ where n cp ~ (T — T c ) d / 2 1 is the number density of fluctuating Cooper pairs. This term is singular for d < 2, but is 
, exactly cancelled by Maki-Thompson (MT) terms. The latter terms are due to new heat transport channels opened 
up by Andreev scattering processes. An electron can Andreev-scatter into a hole, and since electrons and holes carry 
r> \ the same heat current, this leads to a net increase in thermal conductivity. The amplitude for the Andreev scattering 
is exactly the same as for an electron to scatter into a fluctuation Cooper pair, so the MT and DOS terms have the 
same magnitude but opposite sign, and hence cancel. These MT processes lead to a further suppression to electrical 
conductivity since holes carry opposite electric charge to electrons i.e. the MT and DOS contributions cancel for 
thermal conductivity and reinforce for electrical conductivity. 

Before we proceed to the details of our calculation, we present a short history of superconducting fluctuation 
corrections to thermal conductivity. They were first predictedlj in 1970 by Abrahams et al in the diffusive regime. 
These authors concluded that the Aslamazov-Larkin (AL) terms were convergent, but that density-of-states (DOS) 
terms led to divergent contributions in two and one dimensions, of the form In (T - T c ) and (T - T c )- X l 2 respectively. 
They appear to have missed the cancellation between DOSi-apd MT contributions. Shortly afterwards fluctuation 
effects with the predicted power-law behaviour were observedli-3 in one-dimensional Pb-ln wires. After this intial work 
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there was apparently no-theoretical or experimental activity in this area for nearly two decades. Indeed, in Skocpol and 
Tinkham's 1975 reviewo, thermal conductivity is described as one of those quantities which "have not yet benefited 
from sustained interaction between theory and, experiment, perhaps because such effects are small, and hard to 
interpret." In 1990 Varlamov and collaboratorsa predicted AL contributions with the same strong divergence found 
in the electrical conductivity, (T — T c ) d / 2 ~ 2 ; this erroneous result appears to be due to an incorrect treatment of the 
heat-current operator. The same authorso also discussed the relative magnitudes of DOS, MT and AL contributions 
in layered superconductors, and argued that the DOS and MT terms dominate in n c whilst AL terms dominate in 
K Q b. The predicted fluctuation effects have since been seen experimentally in an YE^CuaOy-a single crystaH, and 
Bi 2 Sr 2 CaCu208 and DyBa2Cu307_a polycrystalstffl. Excellent quantitative agreement was found between theory 
and experiment; indeed, even the predicted two- to three-dimensional crossover is seeruat roughly the predicted 
temperature. Fluctuation effects have also been seen in (Nd/Y)BCO intergrowth crystalsB, although these have not 
been compared in detail with theory. However there are problems with this apparent agreement between theocjj and 
experiment. The AL-Gpntributions have been re-analysed in two works using phenomenologieal hydrodynamidl3 and 
Gaussian fluctuation^!! approaches, and argued to be convergent. Very recently Savona et afca have agreed that there 
is no divergent AL correction, but argue that there are still divergent DOS and MT terms; we believe that these 
authors have missed the cancellation between the DOS and MT terms. 

We now proceed to the details of our microscopic calculation. The thermal conductivity is obtained from the 
imaginary time heat response kernel, Qhh(iQn), by analytic continuation from positive Bose Matsubara frequencies, 
fi n = 2itTn, 



k = lim 

f2->o 
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(5) 



The diagrammatic contributions to the heat response kernel of lowest order in perturbation theory are detailed in 
Fig. 1. The solid lines are disordered electron Green functions 
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iei -£k + ^rsgn(e;) 



(6) 



where e% = 2irT(l + 1/2) is a Fermi Matsubara frequency, £/. = fc 2 /2m — /i is the electronic excitation spectrum, and 
r is the elastic scattering time. The black dots represent heat-current vertices, which are given by 



jfc(k, £;,£/+ Q n ) = —i(2ei + fl n ). 
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The shaded regions are impurity vertex renormalization which, at zero momentum, take the form 

e(-eiea) 



C(g = 0,ei,e 2 ) = 0(+£i£ 2 ) + 



(7) 
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whilst the dashed lines are single impurity renormalizations. The wavy lines are superconducting fluctuation propa- 
gators, L(q,iui m ), which for small q are given by 



L(q, iu>r, 
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where N(0) is the electronic density-of-states per spin at the Fermi surface, ui m = ZirTm is a Bose Matsubara 
frequency, ip(x) is the digamma function, D = v^r/d is the diffusion constant, and A(u> m ) is given by 
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The zero-frequency fluctuation propagator, L(g,0), has a 1/g 2 divergence as T approaches T c , 

~T-T C 



L{q,Q)- 1 = JV(0) 



A{0)Dq 2 



(11) 



It is this feature which leads to divergent contributions to various physical properties as T approaches T c . 

Diagrams A and B of Fig. 1, in which a fluctuation propagator affects only one electron line yields the DOS 
contributions; diagram C, in which a fluctuation propagator leads to interference between electron lines, yields the 
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MT contribution; diagram D, which possesses two fluctuation propagators, yields the AL contributions. Note that 
since the object of this paper is merely to show that there are no divergent contributions to k at T = T c , we have 
omitted all diagrams that cannot have such divergences. In particular we have ignored all DOS and MT diagrams 
that have an impurity line or ladder between the two heat current vertices. Such diagrams possess an extra factor of 
q 2 which removes the low-momentum singularity of the fluctuation propagator, L{q, 0). We need also consider only 
the lowest power of q in any diagram since this will have the most divergent behavior - we therefore set q = in all 
terms except the fluctuation propagators. Finally since all DOS and MT diagrams have only one superconducting 
fluctuation propagator we can take the static limit and consider only terms L(q, iuj m ) with zero Cooper pair frequency, 
Lj m = 0. The AL term has two fluctuation propagators, and here we have to be more careful and keep all u> m terms as 
there is an anomalous region of frequencies where one propagator can have positive frequency, and the other negative 
frequency. 

The regular parts of the DOS and MT diagrams, which come from diagrams A and C, give the total contribution 
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The sum of three terms in the curly brackets is easily seen to be proportional to fl^ so upon analytical continuation, 
division by f2, and setting to zero, we get zero contribution. The two DOS and one MT term have exactly cancelled 
each other. Note that the same terms in the electromagnetic response function reinforce rather than cancel each other 
because the electric current vertex has the opposite electron-hole parity to the heat current vertex (ie. holes carry 
opposite charge but the same excitation energy to electrons). 

The anomalous parts of the DOS and MT diagrams give total contribution 
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where we have explicitly carried out the ei sum. Upon analytically continuing ifl n — > f2, dividing by f2, and taking 
the limit SI — ^ 0, the above expression gives zero result. The net result is thus that the anomalous part of the DOS 
+ MT diagrams do not yield a divergent contribution. 

Finally it only remains to show that there is no divergent contribution from the AL terms. Paradoxically, although 
this result does not appear to be in dispute, it is the trickiest to prove. The method used is simple power-counting, 
applied to the analytical continuation of the complete Matsubara sum. We need the complete sum because there 
is an anomalous region of Bose frequency, u> m , for which the two superconducting propagators L(g,«w m + i£l n ) and 
L(q,iuj m ) have opposite signs of Matsubara frequency. We cannot therefore simply take the static approximation 
where one or other superconducting propagator has zero Matsubara frequency. Instead we must evaluate the two 
triangle blocks for general u> rn , and distinguish between the three summation regions: (i) oj m + Q n > 0, u> m > 0; (ii) 
Lo m + Q n > 0, uj m < 0; (iii) w m + Q n < 0, uj m < 0. Note that the two summation terms, w m = and uo m = — Q n , 
which possess one divergent fluctuation propagator, L(q,0), are both zero after analytic continuation iQ n Q + iO, 
division by f2, and taking the limit £1 — > 0. It follows that when we analytically continue using the contours shown in 
Fig. 2., we need not worry about contours passing through the poles. 

The contributions from regions (i) and (iii) give identical results, and their sum is 
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Bi(iu) m , itt n ) 2 L(q, iu: m )L(q, iw m + iVL n ), 
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where the Bi(iuj m , zf2 n ) are from the triangle blocks. Upon replacing summation over u m by integration over u>, and 
analytically continuing iQ n — > S7 + «0, we get 



1 f +oc a 2 
Qf L (n) = -— / dujcoth{uj/2T)^2^-B 1 (uj,n) 2 L(q,uj)L(q,uj + n) 



(15) 



For small u, £1, we can show that Bi(u>, fi) w auj + /3f2, where a and f3 are constants, so that for power-counting 
purposes Eq. (113) at T = T c becomes (ignoring all irrelevant coefficients) 



+ OC 



dwcoth(w/2T) / d d qq 2 



(w + nf 



(q 2 — iui){q 2 — iuj — iff) 



(16) 



3 



The 0(fl) piece can be found by expanding either the numerator or denominator. In both cases the behavior as 
u q 2 ~ is 0(q d ), and hence there is no infrared singularity for d > 0. 
The contribution from region (ii) has the form 

q 2 

Q2 L (itt n ) = -T ^ X! —^ B 2{iu m ,iVL n ) 2 L(q, iuj m )L(q, iu> m + ifl n ) (17) 

0>w m >-£l„ q 



which, upon replacing summation over io m by integration over lo, gives 
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Shifting variable in the second integral, u> — > w — if2 n , analytically continuing if2 n — * il + iO, shifting the variable 
back, u> — > to + Q, dividing throughout by Q, and letting — > gives 

n->o S2 87riT J_ 00 sinh (u)/2T) ^ a 

For small w we can show that i?2(w,0) = 70;, where 7 is a constant, so that for power-counting purposes Eq. (|l^ ) at 
T = T c becomes 

lim ^— i- i - / ? / d d qq 2 — > — 7-5 — 20 

O-rf) J- 00 smh{uj/2T) 2 J (q 2 - iuj)(q 2 + iuj) 

The behavior as w ~ g 2 ~ is 0(q d ), and hence there is no infrared singularity for d > 0. We have therefore shown 
that there is no singular contribution from the AL diagrams. 

In conclusion we have shown that there are no superconducting fluctuation corrections to the thermal conductivity 
above the transition temperature which are singular as T approaches T c . The experimental features seen near T c must 
therefore have some other physical explanation, such as reduced phonon scattering from normal state electrons. We 
hope that there will be continued experimental interest in thermal conductivity near T c in one- and two-dimensional 
superconductors, of both the high-T c and low T c variety. In future work we also intend to evaluate the non-singular 
fluctuation contributions to the thermal conductivity to see if this can explain any of the experimental features 
(although, given their power-law behavior, this seems unlikely). 
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FIG. 2. Contour required to perform sum over Matsubara frequencies u> m in the AL diagram. The branch cuts at 7m(w) = 
and Im(u)) = — Q n come from the fluctuation propagators L(q,iio m ) and L(q,iui m +ifl m ). The poles summed over fall into 
three regions separated by the two branch cuts: (i) uj m > 0; (ii) > uj m > —fin; (iii) —fin > u>m- These contours can be 
deformed to contours parallel to the real axis as shown in the figure. Note that the poles which lie on the branch cuts yield no 
singular contribution and can be ignored. 
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